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Abstract 

The production and decay of new possible heavy Majorana neu- 
trinos are analyzed in hadronic collisions. New bounds on the mixing 
of these particles with standard neutrinos are estimated according to 
a fundamental representation suggested by grand unified models. A 
clear signature for these Majorana neutrinos is given by same-sign 
dileptons plus a charged weak vector boson in the final state. We dis- 
cuss the experimental possibilities for the future Large Hadron Col- 
lider (LHC) at CERN. 



PACS 12.60.-i, 14.60.St 



*marroqui@if.ufrj .br 
' yar a@ if . uf r j . br 
*simoes@.if.ufrj .br 
§aline@if. ufrj.br 



1 



1 Introduction 



The increasing experimental evidence for neutrino oscillations seems to in- 
dicate that physics beyond the standard model may be at sight. The com- 
bined results from solar, atmospheric and laboratory neutrinos, if confirmed, 
require at least one new light neutrino JT|. Besides this possible new fourth 
neutrino, theoretical models must give some explanation on the mass spec- 
trum of neutrinos, as well as on the pattern of mixing angles. This is usually 
accomplished by variations of the see-saw mechanism, with new heavy neu- 
tral leptons. Many models that have being proposed recently could account 
for these properties but all of them need further experimental tests ||. One 
particular point that still needs confirmation is the nature of the neutrino 
fields - Major ana or Dirac. 

In this paper we analyze the production and decay of new heavy Majo- 
rana neutrinos in hadronic collisions. In order to understand the smallness of 
neutrino masses, the most likely mechanism include at least one new heavy 
Majorana neutrino. These neutrinos could have its origin in many grand 
unified models such as in the fundamental 27-plet of Eq or in left-right sym- 
metric models. Here we consider the case of new heavy neutral singlets. In 
practically all these models we have new interactions, new fermions and a 
large Higgs sector. As we have presently no evidence for new interactions, we 
will make the hypothesis that the new neutrino interactions are dominated 
by the standard Z and W weak gauge bosons. New Z' and W are possible 
but they are known to have a small coupling with the standard fermions. 
Since these Majorana fermions can decay in both charged positive and nega- 
tive leptons, they will give us very clear signatures. One interesting signature 
is given by the production of a pair of same-sign dileptons accompanied by a 
charged weak vector boson. This signature was already suggested in models 
with a new right weak boson [0] and composite Majorana neutrinos [|J. In 
a previous work || we have presented estimates for this process in a model 
with heavy Majorana neutrinos mixed with standard neutrinos. Here we 
generalize our results for the mixing of light and heavy neutrinos and present 
new bounds on mixing angles. As we will discuss in section 2, the mixing 
with new possible heavy singlets makes pair production of Majorana neu- 
trinos suppressed relative to a single heavy Majorana production. We have 
calculated the general elementary process q^qj — > £~£~~W + , with £ = e,/z, 
including finite-width effects, and discuss experimental cuts for the LHC at 
CERN energies in section 3 and give our conclusions in section 4. 
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2 The model and bounds on mixing angles 



The key element in our analysis is the mixing between the standard light 
neutrinos and a new heavy one. Recent experimental searches at LEPII 
at CERN imply that these new possible neutral states must have a mass 
above the standard Z mass ||. This kind of mixing is different from the 
light-to-light mixing, which seems necessary in order to explain the standard 
neutrino properties in solar and atmospheric phenomena. We work in a 
scenario where each fermionic family is enlarged and the mixing occurs with 
neutrinos of the same generation. In the naive see-saw model, the mixing 

TTl 

between light and heavy neutrinos is given by 9 ~ and it is clearly out of 

m N 

any visible effect. However, there are many theoretical models which decouple 
the mixing from the mass relation 0. The mechanism is very simple. In 
the general mass matrix including Dirac and Majorana fields one imposes 
some internal symmetry that makes the matrix singular. Then the mixing 
parameter has an arbitrary value, bounded only by its phenomenological 
consequences. We take this mixing to be described by a single parameter 
(for each generation), with the new heavy states appearing as isosinglets. 
The general Lagrangian of this kind of model is then given by: 

C = £ {--^=(^7 M (l-7 5 )/cos^ + iV^(l-7 5 )/sin^)iy / , 

9 /-«/_ _ _.5N.. , 1 



(iVfigv, - 9A l l h )v i + - cos0 t sin fyiV^l - 7 5 )^ 



2 cos 9\y ^ ' 1 ' 2 

sin 2 

2 



\ sin 2 0,%"(1 ->y 5 )N l )z ii + h.c.}, (2.1) 



where 



g Vl = g s v M - 1 -sm 2 9 l 

9 Al = g s A M - 1 -sm*9 l (2.2) 

We have considered only the first and second families, since electrons and 
muons are directly observable particles. The r ± r ± dileptons are also possible, 
but the associated Majorana neutrino can be heavier than the corresponding 
to the other families. Besides that, the final state for tau leptonic decays has 
to be reconstructed with undetected neutrinos. 
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The W and Z couplings to standard leptons are well known to agree with 
the experimental data on low energy lepton-hadron scattering, with the high 
precision Z data at CERN and SLAC, as well as with all the charged cur- 
rent experiments. This imply that the mixing angle Qi must be small. It 
is also well known that the present experimental data requires the standard 
model quantum corrections in order to compare theory and data. With this 
in consideration, we take the changes in the physical observables due to new 
neutrino mixing to be small contributions to the standard model theoretical 
results, including first order corrections. With a em and Mz as fundamental 
input parameters, the mixing indicated in Eq. (2.1) above will imply correc- 
tions to the effective //-decay constant G M ; charged coupling lepton universal- 
ity, Cabbibo-Kobayashi-Maskawa unitarity and neutral current interactions. 
The most stringent bounds on 9 e come from the effective coupling of the Z 
to the electron neutrino § g u e e xp = 0.528 ± 0.085 and Tf^Z) = 500.1 ± 1.8 
MeV, to be compared with the standard model predictions Qsm — 0.5042 
and r^-(Z) = 501.7 ± 0.2 MeV. For the muon neutrino coupling with the Z 
boson, the Particle Data Group quotes = 0.502 ± 0.017. After a global 
fit to the data we have obtained, at 95% confidence level, the upper bounds: 



These bounds are more restrictive than earlier analysis ||. Since the 
muon bound is much stronger, the dimuon channel will be much smaller 
then the dielectron channel for the LHC at CERN. We are left then with the 
possibility of detecting the Majorana contribution in the e ± e ± channel. 

With these bounds we can compare the total cross section for single and 
heavy pair neutrino production. The interaction shown in Eq. (2.1) implies 
that heavy pair production is suppressed if compared with single production 
of a heavy neutrino, as it was shown in 

3 Results 

The fundamental production mechanism is given by quark-antiquark anni- 
hilation into a W followed by the single Majorana production and decay, as 



sin 2 # e < 0.0052 



(2.3) 



sin 2 < 0.0001 



(2.4) 
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shown in Fig. 1. At the LHC at CERN center of mass energies (y/s = 14 
TeV) other contributions such as gluon fusion will give a small contribution. 
We have also checked the contribution from W W fusion which turns out to be 
small. Let us now discuss the main characteristics of the final state identifi- 
cation. Experimentaly each charged lepton can be isolated into a cone, with 
no accompanying hadrons, in order to eliminate possible misidentification 
from hadron decays. The dominant final state W — > jets has a sharp kine- 
matical identification in the invariant mass. The leptonic W — > ivi decay 
can be identified by the total px balance with a reconstructed neutrino. The 
proposed detectors for the LHC at CERN energies are planned to have a 
good leptonic and hadronic resolution, as well as a wide rapidity coverage. 
We have estimated the total cross section for pp — ► £ £ W^X applying 
conservative cuts on the final state. For the charged leptons we considered 
Pti > 5 GeV and a pseudo rapidity \r)\ < 2.5. For the charged weak boson 
the cuts \r]\ < 2.5 and prw > 15 GeV were done. For the upper bound 
sin 2 9 e = 0.0052 our results are shown in Fig. 2. The analytic expression for 
the elementary process — > £~£~W + is given in the appendix. For the 
proton structure functions we have employed the MRS[G] set jlOfl . 

The standard model background comes mainly from top production and 
decay via the chain t — > b — > c £vg. It is well known [|TTJ that prompt lepton 
background from heavy quark decay can be suppressed by a suitable isolation 
criteria on transverse energy and momentum. This procedure was applied 
on earlier analysis [O] of Majorana neutrino production. We employed a cut 
of pt = 80 GeV on the softer lepton coming from the standard model back- 
ground. Since our estimate on mixing angles is more restrictive than previous 
values, the signal for Majorana neutrino production will be suppressed. An 
alternative to enhance the signal is to employ a very distinctive characteris- 
tic between the signal and the standard model background: whereas in the 
Majorana neutrino production we have no undetected particles in the final 
state, for the heavy quark background we have two undetected neutrinos in 
the final state. As a consequence there is no missing transverse momentum, 



P't, for the signal and for the background there is an average |12j p'x — 50 



GeV. Recent analysis [13] on detector possibilities at LHC energies indicate 
that an error on the missing transverse energy in the 10 — 20 GeV region is 
expected. The other kinematical difference between signal and background 
is given by the total hadronic invariant mass in the final state. In the heavy 
lepton production we have the associated production of a single W boson, 
with a corresponding peak at the W mass. The standard model background 
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has a very different hadronic mass distribution. 

We show in Fig. 3 the convolution for the hadronic invariant mass versus 
final state missing energy variables. We have applied a cut on the transverse 
momentum of the least energetic lepton of pt > 35 GeV. Both the signal 
and background events were hadronised according to Pythia fl4j| . The box 
in Fig. 3 shows the kinematical region of the expected signal events. We 
have considered a heavy Majorana neutrino with mass of 200 GeV. From 
Fig. 2 for a total cross section of 1 fb and for the expected luminosity of 100 
fb -1 at LHC we have a total number of 100 events for the signal. For an 
invariant hadronic mass in the 60-85 GeV region, a total missing energy less 
than 12.5 GeV and the most energetic charged lepton transverse momentum 
greater than 35 GeV, we still have 62 events from the signal, whereas for the 
standard model background we have only 8 remaining events. With these 
cuts, we estimate the ratio s/y/s + b, where V and 'b' are the signal and 
background number of events, to be 7.4. This clearly shows the advantages 
of applying the above procedure for the cuts. For higher neutrino masses the 
total cross section fall requires more stringent cuts on the suggested variables. 
In this case other kinematical restrictions may be useful. The standard model 
background can still be reduced one order of magnitude by requiring only two 
hadronic jets in the final state and by increasing the transverse momentum 
cut on the charged leptons. In Fig. 4 we display the transverse momentum of 
the final leptons for some values of the Majorana neutrino mass. For masses, 
in the 100 — 200 GeV region, a high pt cut of 80 GeV strongly reduces the 
signal but for higher masses there are practically no changes. For example, 
for mjy = 200 GeV the cross section shown in Fig. 2 (calculated for px = 5 
GeV) is reduced by a factor three if we increase the px cut to 80 GeV. A 
similar procedure for = 400 GeV reduces the signal by 20% and for 
rriN = 800 GeV there is almost no change. From the expected luminosity 
at LHC, we took = 800 GeV as a superior limit in the neutrino mass 
spectrum. For this extreme value of only one event is expected, with 
no further cuts, for a luminosity of 100 fb' 1 . This is a poor result from the 
statistical point of view, but this signal has a very interesting characteristic. 
In Fig. 5 we show the distribution for the angle between the two final leptons 
in the laboratory frame. We can see that for heavier masses region the two 
leptons tend to be emitted back-to-back. We can then safely expect that 
the background can be reduced by the simple requirement of low missing 
transverse energy in the final state, combined with the hadronic invariant 
mass peaked around the W masss. 
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We have also estimated the cross section for the Tevatron upgrade center 
of mass energies (y/s = 2 TeV) at Fermilab with a luminosity of 1000 pb^ 1 . 
For a new possible heavy Majorana neutrino with a mass of 95 GeV we have 
10 events applying the following cuts: for the leptons and the charged weak 
bosons we used \rj\ < 2.0 and E\ > 10 GeV for the leptons, with a mixing 
parameter sin 2 9 e = 0.0052. This number falls to one for a mass value of 140 
GeV doing the same cuts and using the same mixing angle. 

We now turn our attention to the other kinematical characteristics of the 
final state particles. The final W can be identified in the channel W — > jets 
or in the leptonic channel. The distribution for the total invariant mass 
M$ew = (Pe + Pi + Pw) 2 is shown in Fig. 6. A better identification of the 
heavy neutrino mass is given by the invariant mass M^w- Since we have 
two indistinguishable particles in the final state, we must take the variable 
= (Pi + Pw) 2 including all the combinations between the two leptons 
and the W. This distribution presents a well defined peak at the Majorana 
mass, as shown in Fig. 7. The total width for a 100 GeV neutrino is 0.00065 
GeV, and for a 800 GeV neutrino is 0.058 GeV. As mentioned before, another 
variable that can be useful in the identification of the final state particles is 
the angle between the two charged leptons. The distribution for cosOu is 
shown in Fig. 5, calculated in the hadronic center of mass frame. For lower 
M N mass (100 GeV) they have a larger probability of being parallel. This 
is the case also for the Tevatron at Fermilab with « 90 — 140 GeV. For 
higher masses (800 GeV), they show a preference for being emitted back- 
to-back. Finally we show in Fig. 8 the prw distribution for the final W. 
This is an important point in order to eliminate any possible standard model 
background. For lower neutrino mass a cut not greater than 30 GeV must 
be done in order not to loose the signal. For higher masses the pxw cu t can 
be safely higher. 

4 Conclusions 

We have presented an analysis of the possibility of detecting new heavy Majo- 
rana neutrinos through same-sign dileptons and W detection. The agreement 
between the standard model predictions and the present experimental data 
implies that the mixing angle between heavy and light states must be small. 
For the second family, the upper bound for the mixing angle implies that the 
dimuon final state must be suppressed by one order of magnitude relative to 
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the dielectron channel. For the Tevatron upgrade at ,/s = 2 TeV, we found a 
detectable signal in a small region in the heavy neutrino mass. For the next 
LHC the situation is more favorable, allowing an investigation for new heavy 
Majorana neutrinos in the range 100-800 GeV. The signal to background 
ratio can be increased by a combined cut on the invariant hadronic mass and 
low missing energy in the final state instead of a single high pr cut. 
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Appendix 



The elementary cross section for the process q^qj — > £ £ W + 



u, d, s, c), can be written as pq] : 



d&ij = \j{D - CjdLips, (A.l) 

where 

^ = g^sm^Ml, (A 2) 

D = 2 P14P23 (A 3) 

Pr(-p 1 -p 2 ,M w ,T w ) 2 Pr(p 4 +p 5 ,M Ne ,T m y K ' J 



c _ P24P13 + P14P23 + (~Pl4 - P23 - P24 - gig + gig ~ |A§0gl2 

Pr(-pt - p 2 , M w , T w ) 2 Pr{p± + p 5 , Afjv/, r w )Pr(p 3 + p B , M w , r w ) 

(A.4) 



The propagator effects are included via the function 

Pr(P, M k , r fc ) = (P 2 - M 2 ) - t(T k M k ), k = W, Nt, (A.5) 
with £ = e, /1. 

The four-momenta for the quark, antiquark, primary lepton, secondary 
lepton and boson are respectively pi, p 2 , P3, P4, Ph- To simplify our notation 
we have called the product between two four- momenta as: p mn = p m ■ p n . 
Finally Uij is the standard quark mixing matrix. 
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Figure Captions 



1. Feynman diagram for the elementary process qiq^ — > I I W + . 

2. Total cross section versus M Ne for LHC at CERN and Tevatron at 
Fermilab energies. 

3. Invariant hadronic mass versus total missing momentum for the Stan- 
dard Model (crosses) and signal (bullets) events with M Ne =200 GeV. 

4. Normalized final lepton transverse momentum distributions for Mjv e =100, 
200, 400 and 800 GeV for LHC at 14 TeV. 

5. Normalized dilepton angular distribution cos#« for Mjv e =100, 200, 400 
and 800 GeV for LHC at 14 TeV. 

6. Normalized invariant mass distribution M uw for M7v e = 100, 200, 400 
and 800 GeV for LHC at 14 TeV. 

7. Normalized invariant mass distribution M^v for MAr e =100, 200, 400 
and 800 GeV for LHC at 14 TeV. 

8. Normalized W transverse momentum distributions for MAr e = 100, 200, 
400 and 800 GeV for LHC at 14 TeV. 
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Figure 1: Feynman diagram for the elementary process q^j — > £~£~W + . 
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Figure 2: Total cross section versus M^ e for LHC at CERN and Tevatron at 
Fermilab energies. 
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Figure 3: Invariant hadronic mass versus total missing momentum for the 
Standard Model (crosses) and signal (bullets) events with Mn c =200 GeV. 
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Figure 4: Normalized final lepton transverse momentum distributions for 
Mjv e =100, 200, 400 and 800 GeV for LHC at 14 TeV 
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Figure 5: Normalized dilepton angular distribution cos#« for MN e =100, 200, 
400 and 800 GeV for LHC at 14 TeV. 




Figure 6: Normalized invariant mass distribution M^w for MAr e =100, 200, 
400 and 800 GeV for LHC at 14 TeV 
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Figure 7: Normalized invariant mass distribution M#w for Mjv e = 100, 200, 
400 and 800 GeV for LHC at 14 TeV. 




Figure 8: Normalized W transverse momentum distributions for MAr e = 100, 
200, 400 and 800 GeV for LHC at 14 TeV. 
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